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ABSTRACT 

We have developed a detailed stellar evolution code capable of following the simultane- 
ous evolution of both stars in a binary system, together with their orbital properties. 
To demonstrate the capabilities of the code we investigate potential progenitors for 
the Type lib supernova 1993J, which is believed to have been an interacting binary 
system prior to its primary exploding. We use our detailed binary stellar evolution 
code to model this system to determine the possible range of primary and secondary 
masses that could have produced the observed characteristics of this system, with par- 
ticular reference to the secondary. Using th e luminosities and temperatures for both 
stars (as determined bv lMaund et al.|[2004D and the remaining mass of the hydrogen 
envelope of the primary at the time of explosion, we find that if mass transfer is 100 
per cent efficient the observations can be reproduced by a system consisting of a 15 M Q 
primary and a 14 M Q secondary in an orbit with an initial period of 2100 days. With 
a mass transfer efficiency of 50 per cent, a more massive system consisting of a 17 M 
primary and a 16 M Q secondary in an initial orbit of 2360 days is needed. We also 
investigate some of the uncertainties in the evolution, including the effects of tidal 
interaction, convective overshooting and thermohaline mixing. 
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1 INTRODUCTION 

The supernova 1993J, which occurred in M81, was an un- 
usual event. Its spectrum initially showed hydrogen lines like 
a normal Type II supernova, but these subsequently disap- 
peared giving rise to a Type lb-like spectr um, which shows 
broad helium lines and no silicon lines (Ma theson et al.l 
2000, and references therein). The supernova is thus clas- 
sified as a Type lib. It was suggested that in order for this 
transition from a Type II to a Type lb-like lightcurve to 
happen, only a small amount of the star 's hydrogen envelope 
could remain at the time of explosion ()Podsiadlowski et al.l 
1 19931 : IWoosiev et al. 1994). However, this is difficult to ex- 
plain in terms of the evolution of single stars - a star with the 
observed luminosity of the progenitor of SN1993J would be 
expected to retai n a significant amou nt of its hydrogen-rich 
envelope (see e.g. iDavies et al.1l2008l ). It was therefore sug- 
gested that the presence of a bi nary companion was neces - 
sary to explain the observations (|Podsiadlowski et al.ll 19931 1 . 
In this scenario, the red giant progenitor would fill its Roche 
lobe and be stripped of a lmost all its envelope. 

lAldering et alj (|l994T ) noted that photometry of the pro- 
genitor could not be explained by single star models. Ex- 
cesses in the ultra-violet and in the B-band could best be 
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fitted with a spectrum which was the composite of a red star 
together wi th either a blue star or the net spectrum of an OB 
association. IVan Dvk etall (|2002t) suggested that there were 
nearby blue stars which could account for the blue excess in 
observ ations of the progenitor. Subsequently, iMaund et al.l 
(2004) were able to detect the signature of a massive star at 
the location of the supernova, proving that the supernova 
had come from a binary system. SN1993J is therefore the 
first superno va to have been con firmed to come from a bi- 
nary system. IMaund et al.l (|2004T ) give the parameters of the 
primary as log L/Lq = 5.1±0.3 and logT ff/K = 3.63T0.05, 
while they identify the secondary as having the parameters 
logL/L =5.0 ±0.3 and log T cff /K = 4.3 ± 0.1. 

Initial attempts to model the progenitor system of 
SN1993J, along with the supernova light curves that thes e 
models would produce, were made bv lWooslev et all (1 19941 ). 
These authors concluded that the primary had to be in the 
mass range 13-16 Mq and that at the time of explosion the 
helium core mass was 4.0 ± 0.5 Mq, while the hydrogen-rich 
envelope was around 0.20±0.05 Mq based on the bolometric 
light curve obtained from modelling the explosion of their 
evolution models. 

However, at the time of these models nothing was 
known about the properties of the secondary. Now that we 
have measurements of the luminosity and temperature of 
the secondary, it is fruitful to revisit the modelling of this 
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system to see what new cons traints can be added . Partic- 
ularly, the binary models of IWooslev et ail (|l994h looked 
at systems where the initial mass ratio of the components 
(secondary-to-primary mass) was between around 0.6 and 
0.75. It is clear from the location of the secondary that the 
m ass ratio was much closer to unity. The model presented 
bv lMaund et al.l (|2004T l was an improvement on earlier work 
as these authors had obtained details of the secondary. How- 
ever, the degeneracy of the solution was not studied. Here we 
build upon their model by examining how sensitive the sys- 
tem is to the initial conditions and the treatment of detailed 
physics such as tides, thermohaline mixing and convective 
overshooting. 



2 THE STELLAR EVOLUTION CODE 

We use an extensively modified version of the stel lar evolu- 
tion c ode STARS, which was originally developed bv lEggletonl 
<jl97lt ) and has bee n updated by ma ny authors. The last 
major update was bv lPols et ail ljl995h and the version used 
in this work makes substantial updates to this. Variants of 
Eggleton's code are unique in that they compute the stellar 
structure and chemical evolution simultaneously, iterating 
on all variables at the same time in order to converge a new 
model. Most codes converge a structure and then use this to 
converge the chemical evolution. While this has been shown 
to be unimportant in most phases of stellar evolution, there 
are some situations where si multaneous solut ion may be nec- 
essarjQ (see the discussion in Stancliffe 2006, for further de- 
tails) . The code now follows the chemical evolution of the 7 
energetically important species 1 H, 3 He, 4 He, 12 C, 14 N, 16 
and 20 Ne alongside the structure variables. The inclusion of 
3 He is necessary for the modelling of low mass stars f ollow- 
ing t he discovery of lEggleton. Dearborn fc Lattanziol (2006, 
I2008T I. who showed that a mean molecular weight inversion 
caused by the burning of 3 He leads to mixing on the red 
giant branch. 

Convectio n is treated within th e code using the mixing 
length theory ( Bo hm-Vitensel Il95ct ) and the models were 
evolved with a mixing length parameter, a = 2.0. This 
value is obtained by calibrating to a solar model. Con- 
vective overshooting was employ ed via the prescription of 
ISchroder. Pols fc Eggletonl |l997) and an overshooting pa- 
rameter of 5 ov =0.12 was usecQ The nucle ar reaction rates 
used i n the code were recen tly updated by Stancliff e et al.l 
i|2005h and IStancliffl (|2005l ). as part of an introduction of 
routines for following the nucleosynthesis of around 40 iso- 
topes from D to 32 S and important iron group elements. 



2.1 Binary evolution 

The most substantial modification to the code is the incor- 
poration of the ideas from Eggleton's binary evolution code 



1 Specifically, if mixing, burning and structural changes proceed 
on similar timescales. This is independent of whether one uses a 
Lagrangian mesh or not. 

2 It should be noted that the oversh ooting parameter was only 
calibrated in the mass range 2-7M llSchroder et al.lll997l) . We 
are making the assumption that this value also applies outside 
this mass range. 



twin, (|Eggleton fc Kiseleva-E gglcton 20~02h. The twin code 
simultaneously solves the equations for the evolution of both 
stars in a binary system along with their orbital properties. 
We have modified the STARS code to be able to do this. This 
necessarily requires the inclusion of the physics of binary in- 
teraction. We assume that the orbit of the binary is circular. 
The extension to elliptical orbits is reserved for future work. 
For the systems we are interested in the orbit will circularise 
before the primary fills its Roche lobe. 

In order to model the orbits of the stars, we have added 
three new variables to the code: the angular momentum of 
the orbit and the spin angular momentum of both stars. 
We assume that the stars rotate as solid bodies. Recent 
simulations of stars with both rotation and magnetic fields 
suggest that this may be a reaso nable assumption (at leas t 
for main sequence evolution, see iMaeder fc Mevnetl I2005T ) . 
We account for t idal i nteraction of the binary using the 
equations of iHutl lll98ff) an d we follow the prescriptions of 
iHurlev. Tout fc Pols! (|2002f > for the details of the equilibrium 
tide with convective damping and the dynamical tide with 
radiative damping. 

The code also accounts for the loss of angular momen- 
tum via winds and for the transfer of angular momentum 
from one star to the other, both with direct mass transfer 
from Roche lobe overflow and for wind accretion. We ac- 
count for the variation in composition of the material that 
is transferred from one star to the other. Again, this is done 
both for the case of Roche lobe overflow and wind accretion. 

We find it necessary to modify the physics of accretion, 
mass loss and angular momentum loss in the case that a 
star rotates close to break-up. We do this by modifying the 
mass-loss: 



M(fi) 



M(0) 
n 

ficrit 



(1) 



where M{Q) is the mass-loss rate as a function of rotation 
rate, M (0) is the mass-loss rate in the absence of rotation, SI 
is the angular frequency of rotation and f2 C rit is the angular 
frequency for critical rotation. This increase in M reflects 
the fact that as the star spins more rapidly, its surface layers 
are less gravitationally bound. Also, we presume that a star 
which is rotating at close to its critical rotation rate will 
not continue to accrete more material. We therefore reduce 
the accretion rate by a factor of 1 — n n . Finally, we also 
enhance the rate of loss of spin angular momentum by a 
factor of (1 - ttM" 1 . 

During mass transfer, it is not guaranteed that the ac- 
creting star will be able to accept all the matter than the 
primary transfers. We therefore set a maximum rate of accre- 
tion for a star equal to the star's mass divided by its K elvin- 
Helmholtz timescale, as is done bv lHurlev et~aT1 <|2002h ■ This 
is in addition to the factor discussed above. 

The code also includes the physics for therm o haline 
mixing, using the pr escription of Kippcnh ahn et ail (|l980h 
as implemented by Stancliffe et al. (|2007h . Thermohaline 
mixing is an important process for systems undergoing mass 
transfer. It occurs when the mean molecular weight of the 
stellar gas increases towards the surface. Although this situ- 
ation is dynamically stable as long as the density decreases 
outward, a secular instability sets in when a gas element is 
displaced downwards and compressed. Such an element must 
be hotter than the surrounding gas by virtue of its higher 
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molecular weight and if it loses heat to its surroundings its 
density increases and it will continue to sink. As a result mix- 
ing occurs on a thermal tim escale, until the molecular weight 
difference has disappeared (|Kippenhahn et al]|l980l ). In in- 
teracting binary systems, mass transfer (by wind accretion 
or Roche lobe overflow) may lead to the transfer of mate- 
rial that has been processed by nuclear burning from the 
primary to the surface of the less evolved secondary. This 
alters the surface composition of the secondary. The extent 
to which the surface composition is changed depends on how 
much the transferred matter is mixed with the pristine stel- 
lar material. 



2.2 Opacities 

The opaci ty routines of the code w ere substantially over- 
hauled by lEldridge fc Toutl (1200431 ). These authors modi- 
fied the routines so that they now account for changes in 
the opacity related to variations in the carbon and oxy- 
gen composition of the material. The high temp erature 
opacities are taken from llglesias fc Rogers! i|l996l ). while 
at temperatures below log^ T/ K = 4.0, the opacitie s of 
lAlexander fc Ferguson! (|l994T ) and lFerguson et all (|2005l ) are 
used. Opacity tables from Z = up to Z = 0.05 have 
been produced. The lo w temperature opacities used by 
lEldridge fc Tout! (|2004ah do not account for molecular opac- 
ities. These are only important at low temperatures and are 
particularly important for carbon-rich matter. These condi- 
tions are not encountered in the mass ranges studied in this 
work0 



2.3 Mass-loss rates 

The new version of the code now has some of the most 
commonly used mass los s rates includ ed in it. In addi- 
tion to the Reimers' rate (lReimerslll975h. the c o de als o in- 
cludes the pre scriptions of IVassiliadis fc Wood! l|l993l ) and 
Blocker (1995). These r ates h ave r ecently been employed 
by IStancliffe fc Jeffervl (|2007h and lLau. Stancliffe fc Tout! 
(2008) in modelling the evolution of AGB stars. 

For the massive stars consider ed in this work, w e employ 
the mass-loss rates described in lEldridge fc Toutl (|2004bl ) 
an d lEldridge fc Vinkl d2006h These authors use the rates 
of IVink. de Koter fc Lamersl i|200ll ) for OB stars, which 

, the rates of 



have a metallic i ty sca ling factor of (Z/Zq 
Ide Jager et al.l (Il988 j ) for pre-Wolf-Rayet phases of evolu- 



tion and the rates of lNugis fc Lamersl ((2000) for Wolf-Rayet 
phases. A common metallicity scaling factor of (Z/Zq) ' 5 is 
used for both these phases. 



3 MODELLING AND RESULTS 

Based on the luminosity limits of the primary and secondary, 
we can immediately constrain what mass of star the pri- 
mary was and what the secondary mass would have been 
at the point the primary exploded. If the primary is less 



3 However, Stancliffe fc Glebbeekl | |2008|) have modified the opac- 
ity routines to include the prescription of jMarigo ( 20Q2) i n order 
to study the evolution of asymptotic giant branch (AGB) stars. 



than 11 Mq then the star does not become luminous enough 
to match the lower limit. If it is more massive than about 
32 Mq then the star is too bright to pass through the upper 
luminosity limit. Similarly, the secondary must lie between 
about 14 Mq and 28 Mq . A further constraint is imposed on 
the possible mass ratio of the system. The position of the sec- 
ondary in the HR diagram tells us that it must have reached 
(or be extremely close to) the end of its main-sequence life- 
time and therefore that the mass ratio of the system must 
be close to unity. 

We have constructed a grid of models of varying ini- 
tial stellar masses and initial periods. We use mass pairs of 
12+11 M Q , 15+14 Mq, 17+16 M© and 20+19 M and initial 
periods of 1500, 1750, 2000, 2250 and 2500 days. A metal- 
licity of Z = 0.04 is used as ISmartt et all (120031 ) give the 
metallicity as being roughly twice the solar value. The ini- 
tial abundances use d in the models are solar scaled, based 
on the values from lAnders fc Grevessd (Il989l ). These sys- 
tems are then evolved up to the end of carbon burning. 
The star would explode shortly after this point. We initially 
assume that the mass transfer during Roche lobe overflow 
is 50% efficient. All of the angular momentum lost by the 
mass-losing star is lost along with this material, so that no 
angular momentum is transferred to the accreting star dur- 
ing the mass exchangfl 

The progenitor of the supernova was identified as a 
red supergiant. This means that an interaction between the 
stars had to have occurred when the primary had a convec- 
tive envelope. This presents us with a problem when deal- 
ing with Roche Lobe overflow. Stripping mass from a star 
with a convective envelope leads to the expansion of the star 
and to a positive feedback for mass loss. We find that the 
primary rapidly loses mass and quickly reaches mass-loss 
rates of 10 -2 Mq yr _1 (this val ue is in good agreem ent with 
the mass-loss rates reported bv lMaund et al.ll2004h . Conver- 
gence of the code fails if the mass-loss rate gets much higher. 
We therefore cap the mass-loss rate at this value, purely for 
numerical convenience. It is possible that the mass-loss rate 
in an actual stellar system under these conditions could be 
higher. 

Figure [T] shows the evolution of the 15+14Mq systems 
in the Hertzsprung-Russell (HR) diagram. We find that if 
the initial period of the system is 2250 days or greater the 
secondary is too evolved by the time the primary explodes. 
The model with an initial period of 1500 days encounters 
some numerical problems when the primary star reaches the 
end of core helium burning. The o scillations in the HR dia- 
gram are breathing pulses (see e.g. ICastellani et al ] | 19851 . for 
a discussion of this phenomenon). As core helium burning 
shuts down, the convective core expands slightly, mixing in 
fresh helium. This prolongs the helium burning. It appears 
to be a numerical artefact as increasing the resolution on the 
boundary of the convective core and reducing the timestep 



4 The fraction of the angular momentum that is accreted by the 
mass gaining star presumably depends on the manner in which 
the material is accreted (e.g. through an accretion disc). As we 
do not know what this was in our system we treat it as a free 
parameter. We have run simulations where we allow some angular 
momentum to be accreted by the mass gaining star and they do 
not differ substantially from the simulations presented herein. 
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Figure 1. Example HR diagrams for systems with a primary of 15 Mg and a secondary of 14 Mq. The initial periods of the systems 
are: top left, 1500 days; top right, 1750 days; bottom left, 2000 days; bottom right, 2250 days. Note that the model with the 1500 day 
initial period suffers from breathing pulses at the end of core He burning. The crosses with errorbars denote the observed locations of 
the primary (at logT c ff/K = 3.63) and secondary of SN1993J. 



size at the end of core helium burning reduces the magnitude 
of the effect and in some cases removes it altogether. 

Table Q] gives the final state of the systems evolved. We 
note that the two shortest period systems of the 12+11 Mq 
binary end up with both stars filling their Roche lobes at 
the same time. This system would form some sort of contact 
binary which would require additional physics (e.g. for heat 
transfer between the two stars) that is not included in the 
code. We therefore cannot comment on the outcome of these 
two systems. For the remaining systems, we find that at the 
point of explosion, the primary is barely luminous enough 
to fit within the observed error bars and is also rather cooler 
than the observations suggest. Similarly, the secondaries of 
these systems are also towards the lower end of the errors in 
the observed luminosity. This suggests that the secondaries 
are not massive enough and have not accreted sufficient mass 
during the mass transfer event. The fit could perhaps be 
improved by assuming the mass transfer is more than 50% 
efficient. This is discussed further in section 3.2 below. As 
we increase the initial orbital period from 2000 days up to 
2500 days, the remaining mass in the envelope increases from 
0.66 Mq to 1.90 Mq as the binary widens, with Roche lobe 
overflow setting in later and hence leading to less material 
being stripped. 

The 15+14 Mq systems look more promising. The pri- 



maries all end up close to the observed progentior in terms 
of both temperature and luminosity (as can be seen in Fig- 
ure [1]). The secondaries in the systems with initial periods 
below 2000 days also lie closer to the observed location of 
the secondary, with luminosities lying within the observed 
errors. However, if we have a system with an initial orbital 
period of 2250 days or greater, the secondary is able to evolve 
to a red giant before the primary explodes. This is because 
the wider separation means the primary has to expand more 
in order to fill its Roche lobe. In the time it takes for this 
to happen, the secondary has enough time to evolve off the 
main sequence. Such systems do not reproduced the observa- 
tions of the system. The remaining hydrogen envelope mass 
of the primary varies between 0.07 and 0.65 Mq, with the 
model that has an initial period of 2000 days having a re- 
maining envelope mass of 0.19 Mq. 

The 17+16 Mq systems show the best general agree- 
ment with the location of the progenitor system in the HR 
diagram. The secondaries of all but the system with the 
longest initial period have almost the observed luminosity. 
However, they are all too hot by about 0.2 dex in logT c ff . 
This suggests that the initial mass of the secondary was 
probably slightly higher so that it would have evolved up to 
(or beyond) the end of the main sequence by the time the 
mass transfer took place. The remaining mass in the enve- 
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Table 1. Details of the final state of the models. These models have all been computed with a mass transfer efficiency of 0.5. The 
columns are: the initial masses of the stars in the system in solar masses; initial period of the system in days; the mass of the primary 
at explosion in solar masses; the hydrogen exhausted core mass of the primary at explosion in solar masses; the total remaining mass of 
the hydrogen envelope of the primary at the point of explosion, in solar masses; the logarithm of the luminosity of the primary at the 
point of explosion; the logarithm of the effective temperature of the primary at the point of explosion; the logarithm of the luminosity 
of the secondary at the point of explosion; the logarithm of the effective temperature of the secondary at the point of explosion and the 
N/C ratio in the ejecta material. 



lope of the primary at the time of explosion varies from 0.01 
to 0.44 Mq with the model with an initial period of 2250 
days having 0.20 Mq of envelope remaining at the point of 
explosion. 

For the highest mass pairing, 20+19 Mq, we find that 
mass transfer is too effective at stripping the primary's en- 
velope. Our shortest period systems end up having their 
surfaces completely stripped of hydrogen and so they are 
Wolf-Rayet stars by the time they reach the point of ex- 



plosion. Only the model with an initial period of 2500 days 
comes close to matching the luminosity and temperature 
of the observed progentior. In addition, the secondaries are 
also too luminous and too hot. The former problem could 
be solved by reducing the efficiency of the mass transfer so 
that they accrete less mass. This should make the secondary 
less massive and hence less luminous. 

The constraints of lWooslev et al.l (1 19941 ) allow us to rule 
out several of the models produced above. Their modelling of 
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Figure 2. HR diagram showing evolution tracks for a 14 Mq 
(dotted line), 19 Mq (dashed line) and 28 Mq (solid line) star. 
The observed location of SN1993J's secondary is denoted by the 
cross with error bars. A secondary of around 19 Mq at the time 
of the primary's explosion should give a good match to the ob- 
servations. 
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Figure 3. HR diagrams showing the position of the secondary 
after accretion in a 15+14 Mq system with an initial period of 
2000 days. Each panel displays a track computed with a different 
accretion efficiency, which is marked in the top left corner of the 
panel. The location of SN1993J's secondary is also displayed in 
each panel. 



the SN light curve required a model that has a final hydrogen 
mass of around 0.2 Mq. The 12 Mq progenitors all have final 
envelope masses that are too high and so can be ruled out. 
Of the 15+14 Mq systems, the two shortest period systems 
have too little envelope remaining at the point of explosion, 
while the two longest period systems have too much left. The 
model with an initial orbit of around 2000 days is a good 
match to the Woosley et al. value. The 17+16 Mq system 
with an initial orbit of 2250 days is also a reasonable match. 
The final envelope masses of the 20 Mq systems are all too 
low. 

3.1 Accretion on to the secondary 

In the above models, we assumed that only half of the ma- 
terial transferred from the primary was accreted by the sec- 
ondary. We also noted that in the case of the 12+11 Mq and 
15+14 Mq systems, the secondary is under- luminous com- 
pared to the observations and is therefore presumably less 
massive than is required. We therefore increased the value of 
the accretion efficiency to 1, and re-ran the models of these 
systems. 

In the case of the 12+11 Mq systems, the enhanced ac- 
cretion efficiency meant that the secondary filled its Roche 
lobe during the mass transfer. The code is not equipped to 
deal with these situations and so we terminated the simu- 
lations. Even if we could model this phase of evolution it is 
unlikely that we could form a secondary with a mass suf- 
ficiently large to give it a luminosity of about 10 5 Lq as 
this would require the secondary to have a mass of around 
20 Mq . There is not enough mass in the system to achieve 
this and still have a progenitor with a mass that would fit 
the observations. However, it would be possible to get a sys- 
tem that falls within the error bars of the observation, as 
the lower bound would require only 3 Mq to be accepted 
by the secondary. This is illustrated in Figure [2] Here, we 
plot evolutionary tracks for stars of 14, 19 and 28 Mq. A 
star of less than 14 Mq does not pass through the error bars 
of the observations, nor does one of 28 Mq or greater. We 



note however, that these limits may change slightly when 
the secondary has accreted He-rich material from its com- 
panion as the difference in the opacity of this material will 
shift the star's position in the HR diagram (however, see 
section 13. Tip . 

In the case of the 15+14 Mq runs, the two smallest ini- 
tial periods resulted in the secondaries filling their Roche 
lobes. In the case of the wider period systems, it is possible 
to produce a secondary with approximately the right lumi- 
nosity. However, we note that increasing the mass transfer 
efficiency by a factor of two has little effect on the luminos- 
ity. Despite increasing the final mass of the secondary from 
about 16 Mq up to about 20 Mq , we have only increased the 
luminosity by about 25 per cent. 

The secondaries of the 17+16 Mq systems are of about 
the right luminosity, so we do not vary the accretion ef- 
ficiency for these models. When we look at the 20+19 Mq 
systems, we note that the secondaries are too luminous. If we 
want to reproduce the SN1993J system starting from such a 
high initial mass, we need the accretion efficiency to be lower 
than 0.5 and perhaps the mass transfer has to be completely 
inefficient as a 19 Mq star is of almost the right luminos- 
ity at the point that the primary fills its Roche lobe (see 
Fiugre[2}. However, the simulations show that completely 
inefficient mass transfer does not work: while the primary is 
in the core helium burning phase, the secondary crosses the 
Hertzsprung gap and begins to ascend its red giant branch. 
We would therefore need a slightly less massive secondary 
or a more massive primary in order for this system to work. 

We have also examined the effect of varying the accre- 
tion efficiency for a fixed binary mass and initial separation. 
The evolutionary tracks for the secondary in a 15+14 Mq 
system with an initial orbit of 2000 days are shown in Fig- 
ure [3] If the accretion efficiency is above 0.25 then the sec- 
ondary has a luminosity (and temperature) that are just 
consistent with the error bars of the observations. If the ef- 
ficiency is 1 the secondary has about the right luminosity to 
match the observations, but the star is too blue. We note 
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that if accretion is perfectly inefficient then the secondary 
actually becomes too evolved - accretion is need to slightly 
rejuvenate the secondary. 

We have a similar picture for a 17+16 M Q system with 
an initial period of 2250 days. A complete spread of accre- 
tion efficiencies are consistent with the observed error bars, 
with an accretion efficiency of 0.5 being the best match to 
the observed luminosity. As with the 15+14 Mq system, we 
require there to be some accretion to keep the secondary at 
higher temperatures and prevent it from evolving to the red. 

3.1.1 The effect of thermohaline mixing 

Because the material received by the secondary has been 
processed by nuclear reactions it has a higher mean molec- 
ular weight than the material in the envelope of the sec- 
ondary. We therefore expect thermohaline convection to mix 
this material into the secondary, altering the surface com- 
position. This has not been taken into account in the above 
simulations. 

As a test case, we have re-run the 17+16 M model 
with an initial orbital period of 2000 days. This time, ther- 
mohaline mixing has been included throughout the run. The 
mass transfer is sufficiently rapid that thermohaline mixing 
is unimportant during the mass transfer phase and it is only 
once this phase is over that the effects of thermohaline mix- 
ing become evident. Without thermohaline mixing, the sec- 
ondary has a final surface hydrogen abundance of 0.53 by 
mass fraction, while its helium abundance is 0.42. With the 
inclusion of thermohaline mixing, the final surface H abun- 
dance is 0.61 and the final surface helium abundance is 0.35. 
The final position of the secondary in the HR diagram is not 
greatly affected. With the inclusion of thermohaline mix- 
ing, the secondary is shifted to cooler temperatures by just 
AlogTeff = 0.02, while the luminosity remains unchanged. 
We therefore do not consider thermohaline mixing any fur- 
ther. 



3.2 The effect of tides 

In the above simulations, we have included the effects of tidal 
interaction. As the binaries we are considering are relatively 
wide, we find tidal effects only become important as the 
stars become giants and develop deep convective envelopes. 
This has a particularly important effect on the secondaries 
of many of the systems. During the initially rapid phase 
of mass transfer, the secondary swells up, temporarily be- 
coming a giant. While the mass contained in its convective 
envelope is small, the envelope has a large radius and so the 
convective turnover timescale increases sufficiently to allow 
tidal interaction to influence the system. The spin angular 
momentum of the secondary is increased and as the star con- 
tracts again it spins up, achieving an appreciable fraction of 
its critical rotational velocity. Figure [4] displays this at work 
in a 17+16 Mq binary with an initial period of 1750 days. 

We have constructed a set of models without the inclu- 
sion of tides. The details of the final states of the systems are 
shown in Table [1] Unsurprisingly, the models without tides 
have much higher final masses than the models with tides, 
for a given initial period. Without tides, the two stars do not 
get pulled together as the primary ascends the giant branch 
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Figure 4. Plot of the characteristics of the secondary of a system 
consisting of a 17 Mq primary and 16 Mq secondary during the 
phase of mass transfer. From top to bottom, the panels display: 
the rate of mass gain by the secondary, the mass contained in the 
secondary's convective envelope in solar masses, the logarithm of 
the stellar radius in solar radii and the rotation rate as a fraction 
of the critical rotation rate. The rapid accretion causes the star 
to swell up, temporarily becoming a giant with a convective en- 
velope. Tidal interaction at this point increases the spin angular 
momentum and as the star begins to contract it spins up to an 
appreciable fraction of its critical rotation rate. 

and so a smaller initial period is required to get Roche lobe 
overflow at a given stellar radius. Similarly, the models with- 
out tides have larger masses of hydrogen remaining at the 
point of explosion. The inclusion of tides is equivalent to re- 
ducing the initial period of the system by around 500 days 
for all the mass pairings modelled. For example, in the case 
of a 15+14 Mq binary, a model without tides with an ini- 
tial period of 1500 days produces a primary with similar 
characteristics to a model including tides with an initial pe- 
riod of 2000 days. These cases are shown in Figure [S] In the 
model without tides, the secondary has not had time for it 
to evolve off the main sequence before the primary fills its 
Roche lobe. This results in the secondary being rejuvenated 
by the accretion so it remains blue. 

3.3 Convective overshooting 

All the above models have been evolved using convective 
overshooting as mixing beyond the formal Schwarzchild con- 
vective boundary is needed to reproduce observations such 
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Figure 5. HR diagrams showing the evolution of a 15+14 Mq 
system evolved with and without tides (lower and upper panel 
respectively). The solid line represents the primary and the dot- 
ted line represents the secondary in each panel. The crosses 
with errorbars denote the observed locations of the primary (at 
logT c ff/K = 3.63) and secondary of SN1993J. In the upper panel, 
the system was started with an orbit with an initial period of 1500 
days, while in the lower panel the initial period of the orbit was 
2000 days. 



as the width of the main sequence in clusters (see e.g. 
ISchroder et alj 1 19971 . and references therein) . The mecha- 
nism for this mixing does not have to be actual convective 
overshooting and could be caused by other physical mech- 
anisms such as rotational mixing. We therefore need to ex- 
amine how changing this extra mixing changes the outcome 
of the system. To do this, we have re-run the model grid 
without the use of convective overshooting. 

The details of the final states of the systems are shown 
in Table [21 If we do not include overshooting, we find that 
the 12+11 Mq systems do not undergo Roche lobe overflow 
for any of the periods we have considered. Without over- 
shooting, the stars have smaller core masses and hence are 
less luminous, with smaller radii. While the final remnant 
masses tend to be similar to their cousins with overshoot- 
ing, they have much smaller final He core masses. The non- 
overshooting models lose less of their envelopes, retaining 
larger amounts of H at the point of explosion. 

The models without core overshooting spend less time 
on the main sequence and the end to the main sequence 
is further to the blue. This means that in order to get the 
secondary at the right location it must be a post-main se- 



quence object (as opposed to being close to the end of the 
main sequence) . This places a very serious constraint on the 
possible systems that can produce SN1993J - namely both 
stars must have evolved beyond the main sequence in order 
to reproduce the observations. This is extremely difficult to 
achieve given that the time for crossing the Hertzsprung gap 
is very short. The location of the secondary in the HR dia- 
gram should be considered as evidence for the existence of 
extra mixing on the main sequence. 



4 DISCUSSION 

We find we can adequately reproduce the observed param- 
eters of SN1993J using systems with different initial masses 
and differing accretion efficiencies. A 15+14 Mq system with 
an initial period of 2100 days will produce both a primary 
and a secondary in the right part of the HR diagram if the 
accretion efficiency is 100 per cent. This mo del is in very 
good agreement with the model presented bv lMaund et al.l 
(2004). We obtain a similar primary mass (5.5 Mq compared 
with 5.4 Mq in the Maund et al. model), though our He-core 
mass is lower (4.9 Mq versus 5.1 Mq) and our secondary is 
also of somewhat lower mass (18.8 Mq as opposed to 22 Mq). 
Such a model will also have a final envelope mass of 0.55 Mq 
which is not inconsistent with the observed nature of the su- 
pernova, though this is some what higher tha n the r emaining 
envelope mass suggested by IWooslev et al.l (|l994l ). Higher 
mass systems where the efficiency of mass transfer is re- 
duced can also reproduce the observations. We find that a 
17+16 Mq system in an initial orbit of 2360 days will fit the 
observations provide the mass transfer is 50% efficient. In 
both these cases, the secondary has evolved off the main se- 
quence and would be observed as a blue supergiant. There 
is observational evidence that such sys t ems m ight be quite 
common, as discussed bv lHunter et al.1 l|2008l ). 

The inclusion of tidal physics does not play such an 
important role in determining the final outcome of the sys- 
tems. Regardless of the inclusion of tides, we still need a 
variation in the mass transfer efficiency in order to make 
the secondary of the correct luminosity. 

While comparing the luminosity and surface tempera- 
ture of the models to those dervived from observations is one 
way to determine the quality of our fit, an alternative is to 
calculate a synthetic spectrum for each binary system and 
compare these to the colours observed in the ground based 
pre-explosion ima ges. To do this we use a similar method to 
that described in lEldridge et aLI (120071). We use the BaSeL 
synthetic spectra library ( Westera et al.l ll999T ) to provide 
spectra for each star in our binary. We chose the spectra 
with the surface temperature and gravity closest to that of 
our model star. We then calibrate the spectra so that the 
bolometric luminosity is equal to that in the stellar evolu- 
tion code for a star at 10 parsecs. We work out col ours for 
the U BVRI pass bands using the filter functions of iBesselll 
Jl990l). We set the zero p oint from the flux levels given in 
Drilli ng fc Landoltl (|2000l ) . We hav e taken dust into acc ount 
by using the extinction derived bv lMaund et al.l l|2004l ). 

When we produce the synthetic spectral energy distri- 
butions (SEDs) for our binary systems with this method we 
find that the best fits are the 15 and 14 Mq system with 
a period of 1500 days and the 17 and 16 Mq systems with 
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Table 2. Details of the final state of the models evolved without convective overshooting. These models have all been evolved with a 
mass transfer efficiency of 0.5. The columns are as in Table [TJ 
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periods of 1500 and 2000 days. The most important detail 
to obtain a good match to the observed SED is the surface 
temparture of the red supergiant. However this is the most 
uncertain variable as small changes in the mixing length, 
a, used in the prescription for convection will change the 
radius and surface temperature of the progenitor. Further- 
more such red supergiants may experience oscillations be- 
fore c ore-collapse and th erefore the surface temperature may 
vary (|Heger et al.lll997l ). 

The SEDs of our preferred models do generally match 
the observed SED but not perfectly. For example, for the 
15+14 Mq system which had an initial period of 2100 days, 
we find the red supergiant surface temperature is too low 
with most of the flux coming out redwards of the I-band. 
If we increase logT e ff by 0.15 dex a better fit is found. 
This leads us to favour models where more mass has been 
stripped from the primar y. However, this is con trary to the 
SN lightcurve models of IWooslev et all l|l994h . which re- 
quire a remaining hydrogen envelope of around 0.2 Mq. It 
should also be borne in mind that the theoretical atmo- 
sphere models may not be appropriate to apply to such a 
highly stripped, hydrogen-poor star. 

One may question the extent to which the companion is 
affected by the explosion of the primary star. A substantial 
deposition of energy into the companion and/or stripping 
of material from its surface could shift its position in the 
HR diagram making it hotter and more luminous. The ex- 
tent to which a companion is affected by a SN explosion in a 
binary system has been extensively stu died in potential pro- 
genitor systems fo r Type la supernova ([Marietta et al . 2000; 
iMeng et al.ll2007l ). Our progenitor systems are much more 
widely separated than these Type la progenitors and so we 
would expect the companion to be less affected by the explo- 
sion. There is good evidence to suggest that the companion 
to SN1993J di d not suffer extensiv ely during the explosion. 
The analysis of lMaund et alj (|2004| ) used both pre-explosion 
imaging and post-explosion spectra to determine the nature 
of the companion (see also iMaund fc Smarttl |2009| ) . If the 
companion had been affected by the explosion they would 



have noted a discrepancy. Als o, photometry taken in recent 
years <|Maund fc Smarttl feOOg ) shows that the flux in the U, 
B and V bands (those bands in which the secondary would 
most likely be detected) has retu rned to the same lev els as 
observed prior to the explosion (|Aldering et al.|[l994t ). We 
therefore believe that the companion was not substantially 
affected by the primary's explosion and that its position in 
the HR diagram is a result of its evolution and mass-transfer 
history. 

Additional information would be required to further 
constrain the model. The above assessment is based on three 
observable quantities, namely the location of the primary in 
the HR diagram at explosion, the location of the secondary 
in the HR diagram when the primary exploded and the final 
mass in the envelope of the primary at the point of explo- 
sion. The latter is based on the modelling of the supernova 
light curve and may be regarded as less reliable as it is an 
inferred quantity, rather than a directly observed one. There 
is also a know discrepancy between progenitor masses as de- 
termined from stellar evolu tion calculations and sim ulations 
of Type II SN light c urves (|Utrobin fc Chugaill2008l ). 

iFransson et all (|2005l ) measured the ratio of nitrogen to 
carbon in the ejecta of SN1993J. They obtained a value of 
12.4 + 6. But what does this value represent? If it represents 
the average N/C in the ejecta, then the 15 Mq model falls 
within the error bounds, assuming we eject all the material 
where the binding energy exceeds 10 44 J (i.e. all the material 
that does not collapse to provide the energy for the super- 
nova display). This assumes that the material is completely 
mixed and there is no additional nucleosynthesis in this en- 
velope material during the ejection process. This produces 
4.2 M of ejected material, with an N/C ratio of 10. Ap- 
plying the same criteria to the best 17 Mq model, we get an 
ejected mass of 4.7 Mq and an average N/C of 9.8. However, 
a 20 M primary generates 6.1 Mq of ejecta with an average 
N/C value of 5.4, which is just outside of the observed error 
bar which leads us to favour the lower mass models. The 
reason the more massive star has a lower N/C ratio is that, 
while its H-burning shell is revealed more readily by mass 
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loss, the ejecta contains more C-rich material from below 
this shell. 

We have confirmed that Type lib SN can be the result 
of evolution in a binary system. Whether a system results 
in a Type lib SN is very sensitive to the initial period. This 
suggests that Type lib SN are rare. The range of periods 
over which models produce Type lib SN is around 1000 days. 
If we assume the entire range of initial periods is 3 to 10,000 
days and that the distribution is flat in log 10 P then less 
than 10 percent of binaries give rise to Type lib SNe. 

Recent observations of a second Type lib progenitor for 
SN 2008ax have suggested t hat single stars are als o progeni- 
tors for some Type lib SNe. ICrockett et~al] (|2008f > find that 
this SN likely had a progenitor which was a single Wolf- 
Rayet star with an initial mass of 27 Mq . The mass range of 
such progenitors is small, around 1 Mq . Assuming a Salpeter 
IMF, 1 per cent of single stars will lead to this Type of SN. 
Therefore binaries should be the dominant progenitors of 
Type lib SNe with an upper limit on the relative fraction 
of these binaries being 1 per cent. This figu re agrees with 
the recent observations of lSmartt et al.l (|2008T > who find that 
70.6 per cent of all core-collapse SNe are Type II and of these 
5.4 per cent are Type lib. A detailed prediction of rates is 
beyond the scope of the small set of models presented in 
the paper. However, we suggest that the difficulty found in 
leaving just enough hydrogen on the progenitor leads to the 
objects being rare which is in agreement with observations. 

In producing the model progenitor system we have 
also created a very odd blue supergiant star. It sits in the 
Hertzsprung gap on the HR diagram for 10 percent of the 
main-sequence lifetime and the surface nitrogen abundance 
is enhanced by up to 1 dex. Such stars could e xplain the un- 
expec ted nu mber of blue supergi ants found bv lHunter et al.l 
l|2007h and iHunter et ail (|2008f ). In addition to their en- 
hanced nitrogen abundances, their extended lifetimes as 
supergiants could explain why so many are found in the 
Hertzsprung ga p when few stars shou ld exist there. 

Unlike lPodsiadlowski etail (|l993l ) we find that the com- 
panion star does evolve into a red super giant before explod- 
ing as a SN. iPod siadlow ski et alj (|l993h suggested the com- 
panion of 1993J would explode as a blue supergiant and 
possibly be the progenitor of a SN similar to SN 1987A. We 
suggest this may still be possible if the companion star had 
begun core helium burning before it accreted any mass. In 
our models we find that helium burning has yet to begin and 
the star evolves to become a red supergiant. 



5 CONCLUSIONS 

We conclude that a binary system with a mass pairing of be- 
tween 15-17 and 14-16 M could have given rise to SN1993J. 
The lower en d of this mass r ange i s in good agreement with 
the model of iMaund et al.l |2004). However, there is con- 
siderable degeneracy between the initial mass pairing and 
the efficiency of mass transfer. For an initial mass pairing 
of 15+14 M©, 100 per cent efficient mass transfer is neces- 
sary to produce a secondary in the right location in the HR 
diagram. For a higher initial mass pairing the mass trans- 
fer efficiency can be substantially less, though we believe 
that it cannot be completely inefficient as some rejuvena- 
tion of the secondary is necessary to keep it to the blue. We 



also note that if we were to choose a model based on the 
luminosities and temperatures of the system at explosion 
(without regard to the mass of hydrogen remain as deter- 
mined from lightcurve modelling), we would select models 
th at have smaller hydr ogen envelopes than the predictions 
of IWooslev et~aLl (|l994T ). 

It proves extremely difficult to get the secondary in just 
the right place in the HR diagram. We do not believe that 
this is because the secondary was affected by the super- 
nova explosion. The position suggests that the secondary 
is extremely close to (or just beyond) the end of its main 
sequence. We find that this can only be done in a very nar- 
row range in initial masses and initial periods. Why should 
this binary supernova lie in such a restrictive region? Are 
we missing some important piece of the physics that would 
make it more likely to find companions in this location? 
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